T-dependence of pseudoscalar and scalar correlations 
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Inhomogeneous pseudoscalar and scalar Bethe-Salpeter equations solved using a renormalisation- 
group-improved rainbow-ladder truncation exhibit bound state poles below and above Tc, the critical 
temperature for chiral symmetry restoration. Above Tc the bound state amplitudes are identical, 
as are the positions and residues of the pseudoscalar and scalar poles in the vertices. In the chiral 
limit the it" — > 77 coupling vanishes at Tc, as do /,r, rna-, ga-K-K- For light current-quark masses the 
2ti decay channel of the isoscalar-scalar meson remains open until very near Tc, and the widths of 
the dominant pion decay modes remain significant in the vicinity of the crossover. 
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I. INTRODUCTION 

Numerical simulations of 2-light-flavour lattice-QCD, 
and the study of models that accurately describe dynami- 
cal chiral symmetry breaking and the ir-p mass-difference 
at r = both indicate that a quark-gluon plasma (QGP) 
is reached via a second order transition at a critical tem- 
perature Tc « 0.15 GeV This QGP existed as a 
stage in the early evolution of the universe and its terres- 
trial recreation is a primary goal of current-generation 
ultrarelativistic heavy-ion experiments. The plasma is 
characterised by the free propagation of quarks and glu- 
ons over distances ~ 10-times larger than the proton. 
However, its formation can only be observed indirectly 
by searching for a modification of particle yields and/or 
hadron properties in the debris of the collisions. 

The plasma phase is also characterised by chiral sym- 
metry restoration. Hence, since the properties of the pion 
(mass, decay constant, other vertex residues, etc.) are 
tied to the dynamical breaking of chiral symmetry, an 
elucidation of the T-dependence of these properties is 
important; particularly since a prodigious number of pi- 
ous is produced in heavy ion collisions. Also important is 
understanding the T-dependence of the properties of the 
scalar analogues and chiral partners of the pion in the 
strong interaction spectrum. For example, should the 
mass of a putative light isoscalar-scalar meson [^-Q] fall 
below 2m7r, the strong decay into a two pion final state 
can no longer provide its dominant decay mode. In this 
case electroweak processes will be the only open decay 
channels below Tc and the state will appear as a narrow 
resonance. Analogous statements are true of isovector- 
scalar mesons. 

The Dyson-Schwinger equations (DSEs) Q provide 
a nonperturbative, continuum framework for analysing 
quantum field theories and the class of rainbow-ladder 
truncation DSE models yields a qualitative understand- 
ing of the thermodynamic properties of the QGP phase 
transition at T ^ §| (and also at ^ ||-|ll|). For 
example, using a simple element of this class the pres- 
sure's slow approach to its ultrarelativistic limit, which 



is observed in numerical simulations of lattice-QCD p2[ , 
could be attributed to a persistence into the QGP of non- 
perturbative effects in the quark's vector self energy p^.p] 
Also, using a more sophisticated model, the present in- 
compatibility between lattice estimates of the critical ex- 
ponents [0 was identified as likely an artefact of working 
too far from the chiral limit |l^ . 

The calculation of the T-dependence of hadron prop- 
erties has hitherto used only simple members of this class 
| p^ , p^ 19 1 . Herein we employ a version that provides for 
renormalisability and the correct one-loop renormalisa- 
tion group evolution of scale-dependent matrix elements, 
and focus on scalar and pseudoscalar properties. Mesons 
appear as simple poles in 3-point vertices. Importantly, 
however, these vertices also provide information about 
the persistence of correlations away from the bound state 
pole; e.g., Refs. pO| , ^ , which can be useful in studying 
the T-evolution of a system with deconfinement. Conse- 
quently our starting point is the inhomogeneous Bethe- 
Salpeter equation (BSE). We then employ the homoge- 
neous equation when appropriate and useful. 

In Sec. H we describe the T ^ quark DSE and in- 
homogeneous BSEs, and introduce our model. This also 
serves to make clear our notation. Section III reports our 



results for the T-dependence of scalar and pseudoscalar 
correlations while Sec. IV is a brief recapitulation and 



epilogue. An appendix contains selected formulae. 

II. DYSON-SCHWINGER AND 
BETHE-SALPETER EQUATIONS 

The renormalised quark DSE is 



^T — Q DSE studies characteristically yield momentum- 
dependent vector and scalar quark self energies that remain 
large until fc^ = 1 GeV^; e.g.. Fig. 8 of Ref. @ and Ref J|. 
These features are also observed in lattice simulations um7 
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= 17 ■ P + ^2 («74 i^k + "Ibm) + ^'{PuJt ) : (2) 

where p^^, := (j),ujk) with ujk ~ {2k + 1)ttT the 
fermion Matsubara frequency, and mbm is the Lagrangian 
current-quark bare mass. The regularised self energy is 

(3) 

.A 
Jl,q 

xitr[7'^7^5(9„jr,(g^,;p^J] , (4) 

where: = A,B,C; A,B,C are functions of (|p|^,Ci;^); 

Pa -iZ^/\p\'')n-p.rB ^i,Pc -(Zi/c^fc)i74 ; 

(5) 

and JI^^ :— T J^'il-oo /'^ ^^^/l^"")^; with represent- 
ing a translationally invariant regularisation of the in- 
tegral and A the regularisation mass-scale. The renor- 
malised self energies are 



Hpu>. ; C) = + ^'Ap^. ; A) - ; A) , (6) 

the renormalisation ] 
1 = £,C: and = toh(C). 



C is the renormalisation point, (Cui„)'^ ■= ^ ^-4 = 



A. The Model 

Ti/ (go), ; Pwfc ) in Eq. (^) is the renormalised dressed- 
quark-gluon vertex. It is a connected, irreducible 3-point 
function that should not exhibit light-cone singularities 
in covariant gauges | p2| . A number of Ansdtze with this 
property have been proposed and it has become clear 
that the judicious use of the rainbow truncation 



(7) 



in Landau gauge provides phenomenologically reliable re- 
sults so we employ it herein. A mutually consistent con- 
straint is 



Zi = Z2 and = 



(8) 



The rainbow truncation is the leading term in a l/N^ 
expansion of ; p^^, ) . 

Dfj.^{pn^) is the renormalised dressed-gluon propagator 

g^D^^ipn,) = P/:^{pn,)AF{pnJ + Pj^(poJAG(pn J , 

(9) 

(10) 
(11) 



P^viPn^,) 



P1P3 



fj, and/or = 4, 



Pl^Pu 



Ijl,v = 1, . . . , 4, and Q,k ■= 27rfcT. A Debye mass for the 
gluon appears as a T-dependent contribution to Ap. 

The ultraviolet behaviour of the kernel in Eq. (Q) is 
fixed by perturbative QCD because the DSEs yield per- 
turbation theory in the weak coupling limit. Our model 
is defined by specifying a form for the kernel's infrared 
behaviour: 

AF(PnJ - V{pn,;mg), AG(paJ = 2?(Pn.;0), (12) 
V{pn,;mg) ■.^2Tr^D^dok6^p) + VMipn,;mg), (13) 

where D = (0.881 GeV)^ is a mass-scale parameter and 
47r2 



VMipn^;mg) 



1 - 



In 



(14) 



rrig, T = e — 1, Cij = 1.2TOt, 



with So, = pI^^ 
mt = 0.5 GeV, 7^ = 12/25, = (16/5) tt^T^, and 
Aqcd^ = 0.234 GeV. This model, which is motivated by 
Refs. [^,|^, incorporates the one- loop logarithmic sup- 
pression identified in perturbative calculations. Its pa- 
rameters were fixed at T = by fitting a range of n and 
K meson properties: a renormalisation point invariant 
light current-quark mass rfiu^d = 5.7 MeV, correspond- 
ing to TOfl(lGeV) = 4.8 MeV, gives = 0.14 GeV, 
/tt = 0.092 MeV. The model exhibits a second order chi- 
ral symmetry restoring transition at [ p!?! 

Tc = 0.15 GeV. (15) 



B. Pseudoscalar Channel 

Ward-Takahashi identities relate the 3-point vector 
and axial-vector vertices to the dressed-quark propaga- 
tor. Hence, once a truncation of the kernel in the quark 
DSE has been selected, requiring the preservation of 
these identities constrains the kernel in the BSE. This 
is explored in Ref . pst , where a systematic procedure for 
constructing the kernels is introduced that ensures the 
order-by-order preservation of these identities. 

Using this procedure the inhomogeneous BSE for the 
zeroth Matsubara mode of the isovector ^ vertex con- 
sistent with the rainbow truncation of Eq. (^) is 

r^,(p.,;Fo;C) = ^4iT'75 



Lq 



^g^Dt^APuj^ - qui) 



X l,S{q+ ) Tl^iq^, ; Po; C) S{q-^ ) 7. , (16) 



± _ 



where {t% i = 1,2,3} are the Pauli matrices, g, 

g^^, ± Po/2, Po = {P,0), and Z4 = Zi{C,^) is the mass 
renormalisation constant: 
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mfl(C)r;,(p^,;Pa„;C) 



(17) 



is renormalisation point independent. 

Equation is a dressed-ladder BSE. At T = its 
solution exhibits poles, and their positions and residues 
provide a good description of light vector and flavour- 
nonsinglet pseudoscalar mesons when S is obtained from 
the rainbow quark DSE This truncation is reli- 

able because of cancellations between vertex corrections 
and crossed-box contributions at each higher order in the 
quark-antiquark scattering kernel. 

The solution of Eq. (|l^) has the form (hereafter the 
(■-dependence is often implicit) 



+'i ■ pp ■ P Gl^iPuj,; P) + -fiujk p- PGp^{puj,;P) 



(18) 



where we have neglected terms involving cr^i^-like con- 
tributions, which play a ne glig ible role at T = [ p^ . 
The scalar functions in Eq. (Hq) exhibit a simple pole at 
p2 + = so that 



P2 



■rUPa,,;P)+ regular, (19) 



where "regular" means terms regular at this pole and 
r^(P(j^;P) is the canonically normalised, bound state 
pion Bethe-Salpeter amplitude: 



l,q 



-P) ^^ri(g^,;P)5(,-) 



-f Vl[q^,--P)S{q+)Vl{q^,;P) 



dP 
dP 



, (20) 



P2 = -„12 



with the trace over colour, Dirac and isospin indices, and 
the residue is 



Za tr 



l,q 



(21) 



where X^iQ^ilP) S{q+^)T^{q^^; P) S{q^J is the un- 
amputated Bethe- Salp eter wave function. Substituting 
Eq. (|l^) into Eq. (|6^ and equating pole residues yields 
the homogeneous pion BSE, which provides the simplest 
way to obtain the bound state amplitude. 

At T = 0, r„(C) is the gauge-invariant pseudoscalar 
projection of the pion Bethe-Salpeter wave function at 
the origin in configuration space; i.e, it is a field theoret- 
ical analogue of the "wave function at the origin," which 
describes the decay of bound states in quantum mechan- 
ics. 

The residue of the pion pole in the axial-vector vertex 
is the pion decay constant: 



S^'Pf. 



Z^tr 



l,q 



5^'757xi(gc^,;-P). 



(22) 



It is the gauge-invariant pseudovector projection of 
Bethe-Salpeter wave function at the origin, which com- 
pletely determines the strong interaction contribution to 
the leptonic decay of the pion: 

r^^i., ^i^flGl iKdP ml (1 - ml/mlf , (23) 

Of = 1.166 X 10-5GeV-2, = 0.975, and = 

0.511 MeV, = 0.106 GeV. 

It is a model independent consequence of the axial- 
vector Ward-Takahashi identity that |2^ 



/tt ml = 2m]i{Q r^(C) . 
In the chiral limit 

lim r^(C) = -^{qq)c , 



(24) 



(25) 



where is the chiral limit decay constant, and in this 
model 

/° = 0.088 GeV, -{qqf, goY- = (0.235 GeV)^ 
^ r^(lGeV2) = (0.384 GeV)^ . (26) 



C. Scalar Channel 

The analogue of Eq. ( [l6| ) for the 0++ vertex is pre- 
sented in Eq. (Al). However, the combination of rainbow 
and ladder truncations is not certain to provide a reli- 
able approximation in the scalar sector because here the 
cancellations described above do not occur |2^. This is 
entangled with the phenomenological difhculties encoun- 
tered in understanding the composition of scalar reso- 
nances below 1.4 GeV Q-jl]. For the isoscalar-scalar ver- 
tex the problem is exacerbated by the presence of timelike 
gluon exchange contributions to the kernel, which are the 
analogue of those diagrams expected to generate the -q- 
Vj mass splitting in BSE studies |2^. Nevertheless, in 
the absence of an improved , ph enomenologically effica- 
cious kernel we employ Eq. (Al) in the expectation that 
it will provide some qualitatively reliable insight. (This 
is justified a posteriori.) 



The scalar functions in Eq. (Al) exhibit a simple pole 
at p2 _|_ ^ 



0, Eqs. (A3,A4), with residue 



Za tr 



(27) 



l,q 



where is an obvious analogue of x^- Since a V — A 
current cannot connect a 0^^ state to the vacuum the 
scalar meson does not appear as a pole in the vector 
vertex; i.e, 



ca/3 



PI. 



Z^tv 



2 ' 



t''iX%u>,\P)^Q- 



(28) 



The homogeneous equation for t he sca lar bound state 
amplitude is obtained from Eqs. (Al,A3). 
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D. Two-body Decays 

The a and tt bound state amphtudes along with the 
dressed-quark propagators are necessary elements in the 
definition of the impulse approximation to hadronic ma- 
trix elements. For example, the isoscalar-scalar-TTTr cou- 
pling is described by the matrix element in Eq. (A5). 
This yields the width 



16 IT mo 



(29) 



which vanishes if < 2mT^. A contemporary analysis 
of TTTT data identifies a uu+dd scalar with |^ 



0.46 GeV, - 0.22-0.47 GeV, 



(30) 



which corresponds to g^TTTr ^ 2.1-3.0 GeV = 4.5-6.6 mo-. 

We also consider the electromagnetic decay of the neu- 
tral pion, for which the dressed-quark-photon vertex is 
also required in calculating the impulse approximation to 
the coupling. Quantitatively reliable numerical solutions 
of the T = vector vertex equation are now available pO 
However, this anomalous coupling is insensitive to detai 
and an accurate result requires only that the dressed ver- 
tex satisfy the vector Ward-Takahashi identity. This is 
similar to the pion form factor: the = value is fixed 
by current conservation and only the charge radius re- 
sponds to changes in the vertex |Q, and the 7*7r° — > 7 
transition form factor whose value is fixed at the real 
photon point [ pO| but whose g^-evolution is sensitive to 
details of the vertex | |3l| ]. 

An efhcacious T = vertex Ansatz is given in Eq. ( A6) 
p2| , ^ and using this in calculating the m = 7r° ^ 77 
couphng (Eq. (^) for T ^ 0) one obtains 



(31) 



independent of the model parameters. Using this chiral 
limit coupling on-shell 



' ^ 



647r 



(32) 



i.e., 7.7keV cf. the experimental value |33]: 7.7 ±0.6. 



III. MESON PROPERTIES 

Solving the (in) homogeneous BSE at T = is a de- 
manding numerical task because an accurate solution re- 
quires a large amount of computer memory and/or time. 
These problems are exacerbated at T 7^ because of 
the loss of 0(4) invariance, which is manifest in the sep- 
aration of the four-momentum into a three-momentum 
and a Matsubara frequency. Hence in the sum over 
fermion Matsubara modes we usually limit ourselves to 



I = —4, . . . , 4. At the critical temperature this corre- 
sponds to uji^4 > 4.0 GeV, which is greater than the other 
mass-scales in the problem, and yields results that are 
numerically accurate to within ^ 5%. Reproducing the 
T — limit requires many more Matsubara modes 
which precludes that as a check of our numerical method. 
Instead we estimate the error by reducing the number of 
modes and comparing the results. Our discretisation of 
the three-momentum grids in the quark DSE and BSE 
is less seriously limited, and we typically employ > 1000 
Gaussian quadrature points. 



A. Chiral limit 

For TO = t he T — analogues of the inhomogeneous 
BSEs, Eqs. (|l|,^, exhibit poles at 



= and TO^ = 0.56 GeV, 



(33) 



with Too- — 0.59 GeV at to = 5.7 MeV. A low-mass scalar 
is typical of the rainbow-ladder truncation. However, 
there is some model sensitivity; e.g., cf. this result with 
m„ = 0.59 GeV in Ref. to,, = 0.67 GeV using the 
model of Ref. |2£| and Too- = 0.72 GeV in the separable 
model of Ref. [ 37 1 . These four independent calculations 
give an average-TOg. = 0.64 GeV with a standard devia- 
tion of 10%. The rainbow-ladder truncation yields de- 
generate isoscalar and isovector bound states, and ideal 
flavour mixing in the 3-flavour case. Hence the distri- 
bution of mass estimates between that of the isoscalar 
cr and isovector ao(980) might be anticipated. In con- 
trast the three comparison studies, Refs. [ p4p^ , |37| ] , give 
iTiui — mp = 0.75 GeV with a standard deviation of < 2%, 
illustrating the dependability of the truncation in the vec- 
tor channel. 

As already remarked, improvements to the kernel are 
required in the scalar channel. In the isoscalar-scalar 
channel, because Fct/too- is large, it may even be nec- 
essary to include couplings to the dominant tttt mode, 
which can be handled perturbatively in the uj-p sector 
Ipl,^ . In the absence of such corrections, the a properties 
elucidated herein are strictly only those of an idealised 
chiral partner of the tt. Hitherto no model bound-state 
description escapes this caveat. 

The evolution with T of the pole positions in the solu- 
tion of the inhomogeneous BSEs is illustrated in Fig. 
from which it is clear that: 1) at the critical temperature, 
Tc, we have degenerate, massless pseudoscalar and scalar 
bound states; and 2) the bound states persist above Tc, 
becoming increasingly massive with increasing T. These 
features are also observed in numerical simulations of 
lattice-QCD |l|l. 

We obtain the bound state amplitudes from the homo- 
geneous Bethe-Salpeter equations, which Fig. |l| demon- 
strates are certain to have a solution. Their T-evolution 
is depicted in Fig. |^, which indicates that: 1) in both 
cases all but the leading Dirac amplitude vanishes above 
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Tc; and 2) the surviving pseudoscalar amplitude is point- 
wise identical to the surviving scalar one. These re- 
sults indicate that the chiral partners are locally identical 
above Tc, they do not just have the same mass. 



S 60 - 



S 40 - 




FIG. 1. Leading Dirac amplitude for the pseudoscalar 
(shaded sym bols) and scalar (open symbols) vertices, E in 
Eqs. (^,A2), evaluated at {p = 0,iuo) and plotted as function 
of in the chiral limit; i.e., E(p = 0,u!o;P^), for a range 
of temperature values. The bound state poles are evident in 
each case. (We use = 19GeV throughout.) 
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FIG. 2. Dirac amplitudes characterising the m = pseu- 
doscalar and scalar bound states, evaluated at (p = 0,ujo). 
G^r, and Fa- behave similarly: as T — » T~ their value at 
{p — 0,u!o) increases rapidly while the domain on which the 
functions are nonzero shrinks quickly. Esym is the amplitude 
calculated in the chirally symmetric Bo = phase: above Tc, 
Ett = Ecr = Ssym whilc all other amplitudes vanish. 

It is easy to understand this algebraically. The BSE 
is a set of coupled homogeneous equations for the Dirac 
amplitudes. Below Tc each of the equations for the sub- 
leading Dirac amplitudes has an "inhomogeneity" whose 
magnitude is determined by Bq, the scalar piece of the 
quark self energy which is dynamically generated in the 



chiral limit. Bq vanishes above Tc eliminating the in- 
homogeneity and allowing a trivial, identically zero so- 
lution for each of these amplitudes. Additionally, with 
Bq = the kernels in the equations for the dominant 
pseudoscalar and scalar amplitudes are identical, and 
hence so are the solutions. 

It follows from these results that the Goldberger- 
Treiman-like relation ETll 



(34) 



is satisfied for all T only because both /° and Bo{p^^) are 
equivalent order parameters for chiral symmetry restora- 
tion. This possibility was overlooked in Ref. Fur- 
ther, above Tc, the other constraints on the chiral- limit 
pion Bethe-Salpeter amplitude derived in Ref. from 
the axial-vector Ward-Takahashi identity are trivially 
satisfied. 

Using the bound state amplitudes and dressed-quark 
propagators we calculate the matrix elements discussed 
in Sec. H. Their chiral limit T-dependence is depicted 
in Fig. ^ As indicated by Fig. |l|, the pseudoscalar and 
scalar bound states are massive and degenerate above Tc- 

Below Tc the scalar meson residue in the scalar ver- 
tex. To- in Eq. (27), is a little larger than the residue of 
the pseudoscalar meson in the pseudoscalar vertex, in 
Eq. (|2^). However, they are nonzero and equal above 
Tc, which is an algebraic consequence of Bo = and the 
vanishing of the subleading Dirac amplitudes. As a bona 
fide order parameter for chiral symmetry restoration 



A cx (1 - T/Tcf , T/Tc < 1 



(35) 



where (3 is the zero-external-field critical-exponent for 
chiral symmetry restoration (f3 — 1/2 in rainbow-ladder 
models) 0. U^O and ^ for T > Tc demonstrates 
that the pion disappears as a pole in the axial- vector ver- 
tex js^ but persists as a pole in the pseudoscalar vertex. 
Our analysis yields 



oc (1 - T/Tcf , T/Tc < 1 : 



(36) 



within numerical errors, which is evident in Fig. ^j: 
follows a linear trajectory in the vicinity of Tc- Such 
behaviour in the isoscalar-scalar channel might be antic- 
ipated because this channel has vacuum quantum num- 
bers and hence the bound state is a strong interaction 
analogue of the electroweak Higgs boson ||^ . 

m^yj^ in Fig. || is the mass obtained when the chirally 
symmetric solution of the quark DSE is used in the BSE. 
{Bq = is always a solution in the chiral limit.) For 
T > Tc, mly^{T) is the unique meson mass-squared tra- 
jectory. However, for T < Tc, rn^ym < 0; i-e., the solution 
of the BSE in the Wigner-Weyl phase exhibits a tachy- 
onic solution (cf. the Nambu-Goldstone phase masses: 
> toJ = 0) . By analogy with the u-model this tachy- 
onic mass indicates the instability of the Wigner-Weyl 
phase below Tc- It translates into the statement that the 
pressure is not maximal in this phase. 
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0.3 



associated homogeneous equations. Their T-dependence 
is characterised in Fig. ||, which shows that for m 7^ 
they are locaUy identical for T > ^Tc- 
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FIG. 3. m = results. Upper panel: T-dependence of 
the meson masses-squared and pole-residue matrix elements, 
Eqs. (E]JJ22U27^) , along with the mass-squared: rrisy^, and 



residue: 



J, calculated in the chirally symmetric Bq = 



phase. For T < Tc., 



< 0, which is a signal of the insta- 



bility of the chirally symmetric phase at low T. For T > Tc. 

~ ^CT ~ ^sym- LoWel 

and pion decay constant, 



— — mgym- Lower panel: T-dependence of the masses 



within numerical error. 



Figure ^ depicts the evolution of the (common) meson 
mass at large T. As expected in a gas of weakly inter- 
acting quarks and gluons 



mi, 



2 [Jo 



(37) 



where ujq — -kT is a quark's zeroth Matsubara frequency 
and "screening mass." 



B. Nonzero light current-quark masses 

It is st raightforward to repeat the calculations of 
Sec. Ill A for nonzero current-quark masses where chiral 
symmetry restoration with increasing T is exhibited as a 
crossover rather than a phase transition. The solutions 
of the inhomogeneous BSEs again exhibit a pole for all T 
and we determine the bound state amplitudes from the 
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FIG. 4. T-dependence of for large T with m = 0, see Fig. 
= rria for T > Tc and m/{2nT) 1~. This behaviour 
persists with m ^ 0, however, providing an illustration re- 
quires significantly more computer time. (That is easier in a 
simpler model |ig,|l§.) 
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FIG. 5. Dirac amplitudes characterising the m 7^ pseu- 
doscalar and scalar bound states, evaluated at (p = 0,ujo) and 
plotted as a function of temperature. 

Figure H is the rh ^ analogue of Fig. ||. That the tran- 
sition has become a crossover is evident in the behaviour 
of /tt. The meson masses become indistinguishable at 
T ^ 1.2 Tc, a little before the local equivalence is mani- 
fest in Fig. 1^, which is unsurprising given that the mass 
is an integrated quantity. The small difference between 
To- and below Tc is again evident and they assume a 
common value at the same temperature as the masses. 

Figure ^ illustrates the preservation of the axial- 
vector Ward-Takahashi identity via the mass formula of 
Eq. (|2^). The magnitude and T-dependence of both sides 
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are equal within numerical errors above and below the 
crossover. mniC,) r^iC) is the renormalisation point inde- 
pendent quantity that appears in the pion's current alge- 
bra mass formula: nifi^Q rJ(C) and mj^(C) ''7r(C) differ by 
< 5% until T > 0.95 Tc. This comparison illustrates the 
T-domain on which the current algebra formula is valid 
and the analysis of Ref. |^^. It was noted in Ref. | |T7| ] 
that 

r1{C)/f° « (1 - T/T,)-^ , T/T, < 1 , (38) 
which is qualitatively apparent from Figs. ^ and |^. 
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FIG. 6. T-dependence of the meson masses and 
pole-residue matrix elements, Eqs. (|2ll,b2|j27[), for m 7^ 0. 



C. Triangle Diagrams 

We also studied the matrix ele men ts describing the 
two-body decays discussed in Sec. 11 D. With every ele- 



ment in the calculation only known numerically this too 
is a challenging numerical exercise, which we simplified 
by a judicious choice of the external momenta. 

The chiral limit isoscalar-scalar-TTTr coupling and width 
obtained from Eq. (A5) are depicted in Fig. ^, which 
indicates that both vanish at Tc in the chiral limit. Again 



this can be traced to Bq — > 0. For m ^ 0, the couphng 
reflects the crossover. However, that is moot because the 
width vanishes just below Tc where the isoscalar-scalar 
meson mass falls below 2m^ and the phase space factor 
vanishes (see the lower panel of Fig. ||) . 

0.0016 I ' ; ' ; ' ; ' ; ' 



0.0014 




0.0010 



0.14 0.16 
T [GeV] 



FIG. 7. T-dependence of the residue of the pion pole in 
the axial- vector vertex, f-n-m^, and in the pseudoscalar ver- 
tex, 2 m_R((^) r,r(C)- They are equal within numerical errors 
as required by the axial-vector Ward-Takahashi identity. For 
comparison we also plot 2m_R((") r'^{C) and the residue of the 
scalar meson pole in the scalar vertex, 2 m_R((") ro-(C)- 



0.1 




0.14 0.16 
T [GeV] 



FIG. 8. T-dependence of the isoscalar-scalar- tttt cou- 
pling and width, both in the chiral limit and for real- 
istic light current-quark masses. The phase space factor 
(1 - (2m^/m^)^)^/^ is e{Tc - T) in the chiral limit but non- 
trivial for m 7^ 0, vanishing at T ~ 0.98 Tc; i.e., this decay 
channel closes at a temperature just 2% less-than Tc. 

The T-dependence of the 7r'^77 coupling, which sat- 
urates the Abelian anomaly at T = 0, is calculated 
from Eqs. ([A£ , A12 ) and depicted along with the width 
in Fig. p. In the chiral limit the width is identically 
zero because m^r = and the interesting quantity is: 
r(0) = g^o^^lfl- Clear in the figure is that r(0) van- 
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ishes at Tc. It vanishes with a mean field critical expo- 
nent, as is most easily inferred from Fig. |l0[ (An accurate 
calculation is possible because Eq. (U) obviates the need 
for a solution of the BSE.) Thus, in the chiral limit, the 
coupling to the dominant decay channel closes for both 
charged and neutral pions. These features were antici- 
pated in Ref. Further, as is evident in Fig. our 
calculated T(0) is monotonically decreasing with T, sup- 
porting the perturbative 0{T^/ f^) analysis in Ref. H]. 



simulations of lattice-QCD, and find m^^Tr 
T —> oo. 



2'kT as 
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10 
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FIG. 9. T-dependence of the coupling T(0) in Eq. ( |Ai; 
and the vr" — > 77 width. 



For m 7^ both the couphng: g-KO-y-y/ Jtt, and the width 
exhibit the crossover with a slight enhancement in the 
width as T ^ Tc due to the increase in TOtt- There are 
similarities between these results and those of Ref. |Q al- 
though the T-dependence herein is much weaker because 
our pion mass approaches twice the T 7^ free-quark 
screening- mass from below, never reaching it, Eq. (37); 
i.e., the continuum threshold is not crossed. 



IV. SUMMARY AND CONCLUSION 

We employed a renormalisation group improved 
rainbow-ladder truncation of the quark Dyson-Schwinger 
equation, and pseudoscalar and scalar Bethe-Salpeter 
equations to estimate the T-dependence of a range of 
properties that characterise correlations in these chan- 
nels. The rainbow-ladder truncation is quantitatively 
reliable in the pseudoscalar channel. However, that is 
not certain in the scalar channel where the cancellations 
that assist in the pseudoscalar channel are not apparent. 
Nevertheless, we anticipate that many of the features we 
exposed in the scalar channel are qualitatively correct. 

The solutions of the inhomogeneous Bethe-Salpeter 
equation (BSE) exhibit poles at all T, both above and 
below the critical temperature, and in the chiral limit 
and for realistic light current-quark masses. We use the 
associated homogeneous BSEs to determine the masses, 
which correspond to the screening masses determined in 



4 DSE B„ 
8 DSE B„ 
;„,.=16 DSE B„ 
;„„=32 DSE B„ 
;„..= 4 BSE f, E, 




0.13 0.14 
T [GeV] 



0.16 



FIG. 10. T-dependence of gl^^ = /"T(0). Within numer- 
ical error, gl^^ cx (1 - T/Tc). Near Tc, f^(x{l~ T/Tcf^'^ 
because our model has mean field critical exponents and hence 
T(0) (X {l-T/Tcf^. 

The BSE solutions indicate a local equivalence be- 
tween the isovector-scalar and -pseudoscalar correlations 
above the chiral restoration transition/crossover ;0 i.e., 
the scalar functions characterising the bound state am- 
plitudes are identical, and from this follows equality of 
the masses and many of the matrix elements.^ Further, 
the axial-vector Ward-Takahashi identity and the pseu- 
doscalar mass formula that is its corollary are valid both 
above and below the crossover. 

For realistic light current-quark masses the isoscalar- 
scalar meson mass does not fall below 2m^ until very near 
the transition temperature. Hence this dominant decay 
channel remains open almost until the phase boundary 
is crossed. Once it is crossed, however, only electroweak 
decay channels are open. In the chiral limit the anoma- 
lous two-photon coupling constant vanishes at just as 
does /tt, the coupling that determines the strength of the 
leptonic charged pion decay. For realistic masses, how- 
ever, the widths for the leptonic charged pion mode and 
two-photon tt" mode remain significant in the vicinity of 
the crossover. 



^We refrain from asserting this of the isoscalar-scalar cor- 
relation because we cannot anticipate the efi'ect of timelike 
gluon exchange contributions present in the BSE kernel in 
this channel. However, they are kindred to those encountered 
in analysing the realisation of C/a(1) symmetry above the chi- 
ral transition |l|,^,p9| . 

■^The evolution to equality of the masses has been observed 
in numerical simulations of lattice-QCD and in the explo- 
ration of other models that accurately describe dynamical 
chiral symmetry breaking. 
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The construction of a DSE-BSE truncation that al- 
lows for an improved description of the scalar channel 
at r = would provide the foundation for a significant 
improvement of our analysis. Employing an Ansatz for 
the kernel of the quark DSE whose infrared form exhibits 
some T-dependence, perhaps constrained by lattice sim- 
ulations of the string tension |^], may also be interesting. 
However, given the results of Refs. we do not ex- 

pect such a modification to have a significant qualitative 
impact. In the absence of such improvements we never- 
theless expect the local equivalence we have elucidated to 
be exhibited by all isovector chiral partners in the strong 
interaction spectrum. However, the explicit demonstra- 
tion of this is difficult; e.g., in the p-ai complex the bound 
state amplitudes have eight independent amplitudes even 
at r = compared with the four in the pseudoscalar and 
scalar amplitudes at T ^ 0. 



ACKNOWLED GMENTS 

We acknowledge interactions with D. Blaschke and 
Yu.L. Kalinovsky. C.D.R. is grateful for the support 
and hospitality of the Special Centre for the Subatomic 
Structure of Matter at the University of Adelaide dur- 
ing a visit in which some of this work was conducted, 
and both C.D.R. and S.M.S. are grateful for the same 
from the Physics Department at the University of Ros- 
tock during a joint visit in which aspects of this work were 
completed. S.M.S. acknowledges financial support from 
the A.v. Humboldt foundation. This work was supported 
by the US Department of Energy, Nuclear Physics Divi- 
sion, under contract no. W-31-109-ENG-38, the National 
Science Foundation under grant nos. INT-9603385 and 
PHY97-22429, and benefited from the resources of the 
National Energy Research Scientific Computing Center. 



(A2) 



(NB: Here the requirement that the neutral mesons be 
charge conjugation eigenstates shifts the p- P term cf. 
the amplitude.) 



The scalar functions in Eq.(Al) exhibit a simple pole 
at P^ +ml = 0: 



p2 



■r;(p^,;P) + regular, (A3) 



where r^(p(^^;P) is the canonically normalised, 0+^ 
bound state Bethe-Salpeter amplitude: 



2S">^P = tr 



dP 
dP 



(A4) 



The dressed-quark propagator and canonically nor- 
malised Bethe-Salpeter amplitudes make possible the def- 
inition of the impulse approximation to the isoscalar- 
scalar-TTTT matrix element: pi^ — —rn^ — P2^, {p = 



Pi +P2f = -ml, 



9a-^7r ■■= {n{pi)n{p2)\(7{p)) = 
2N,tTD f rAk^,;p)Su{k++) 

Jl,q 



xiT^{kQ+; -pi) S'n(/c+-) ir^(fc_o; -P2) Su{k — ) , (A5) 

kap — k^i + (a/2)pi -I- (/3/2)p2, with only the trace over 
Dirac indices remaining. 



APPENDIX: COLLECTED FORMULAE 



2. Neutral Pion Decay 



1. Scalar Vertex 

The inhomogeneous ladder-like Bethe-Salpeter equa- 
tion for the zeroth Matsubara mode of the 0++ vertex 



ig^D^i^iPu, - qui) 



Lq 



Xl^^S{q+)^liq^,;Pn„■X)S{q-U,, 



(Al) 



where a = 0,1,2,3 with t° = diag(l,l). (NB: In this 
truncation the isoscalar and isovector states are degen- 
erate, which exemplifies our observation that the ladder- 
like truncation is accurate for vector mesons but requires 
improvement before it is quantitatively reliable in the 
0++ sector.) The solution has the form 



An efficacious Ansatz for the dressed-quark-photon 
coupling at T = is 



«r^(^i,4) = iSA(^?,^2)7M (A6) 
+(^i+4)m [^,ii-iei+i2)AAiilel) + ABiilil)] ; 
^F{ilJl)^UF{el) + Fi£l)], (A7) 



(A8) 



where F — A,B; i.e., the scalar functions in the dressed- 
quark propagator: S~-^{p) = 17 ■ pA{p'^) + B{p'^), so that 
this model is completely determined by S(j)). Improve- 
ments of this Ansatz, such as those canvassed in Ref. 
do not have a qualitatively significant effect in the present 
context. 

The dressed-quark-photon vertex is a necessary ele- 
ment in the calculation of the impulse approximation to 
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the vr" — > 77 amplitude. At T = the anomalous con- 
tribution to the divergence of the axial-vector vertex is 
saturated by the pseudoscalar piece of the pion Bethe- 
Salpeter amplitude [Q 

/•A 

T^u{kuk2) = tr / S{qi)^^ThE^{q-~P) 

Jl,q 

X Siq2)tQT^{q2,qi2) S{qi2) tQTA<li2, Ql) , (A9) 

where Q = diag(2/3, —1/3) and herein ki — (fci, 0), k2 — 
(^2,0), P ^ ki + k2, qi = q^ji - fci, q2 = q^^, + k2, q = 
^(gi+gg), 912 = q^j, -fci + fc2- Using Eq. ( [A6| ) to evaluate 



Eq. ( [A9| ) for real photons at T = one obtains 

2)1^(^1,^2) = ^^Cpi/pcr fclp feo- T(0) , 
TT 

with in the chiral limit |30j 

/Or(0) .9,0^^ = 1/2 . 



(AlO) 



(All) 



At nonzero T the tensor structure of Eq. (AlC) survives 
to the extent that, with our choice of ki, ^2, it ensures one 
of the photons is longitudinal (a plasmon) and the other 
transverse. In this case we determine the T-dependence 
using 



Ti4(fcl,fe) 



(fci X k2)^T{0) 



(A12) 



and a generalisation of Eq. (A6) to nonzero T: 

«r(9a;,, , Qc^i., ) = ^Aiql,^ , ql,^ ) n 

+ iqi + q2)[^tG{q^,^,q^J + ABiql^^,qlj], (Al3) 
i^iiq^^t, : ) = ^ciql,^ , ql,^ ) in 

+ {u:i,+uJi,)[\iG{q^,^,q^,^)^AB{ql,^,qlJ], (A14) 
G(g^,^,g^,J =7 - (91 +q2)AA{q\,q\) 

+jiiu;i,+ui,)Ac{ql^^,q^^J, (A15) 

which satisfies the vector Ward-Takahashi identity 

(gc^i, - quji.,)t,iT^{qu,,^,quj,.^) = 5"^(g^,J - S^^iq^,^). 

(A16) 
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